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been  addressed.  The  examples  Investigated  suggest  that  such  shapes  can  be 
found.  Other  results  demonstrate  that  surface  smoothing  of  the  pressure  will 
\  not  be  synonomous  with  the  wave  drag  elimination  due  to  focusing  and  envelope 
formation  of  the  wave  system  off  the  body. 

^For  the  wall  Interference  portion  of  the  effort,  the  method  of  matched 
asymptotic  expansions  Is  utilized  to  study  the  singular  perturbation  problem 
relevant  to  transonic  airfoils  confined  by  large  height  to  chord  ratio 
solid  walU.  In  this  class  of  flows,  the  (Inner)  near  field  Is  represented 
as  a  linear  perturbation  about  the  nonlinear  free  field  which  Is  assumed  to 
be  governed  by  the  Karman-Guderley  small  disturbance  theory  which  Is  non- 
uniformly  valid  as  the  walls  are  approached.  In  the  far  field  (outer)  region, 
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In  which  the  Inner  limit  of  the  outer  solution  is  employed  as  a  far  field 
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incidence  to  obtain  an  Interference- free  value  for  the  lift  are  demonstrated 
from  the  examples.  By  virtue  of  the  nature  of  the  perturbation  method,  the 
height  dependence  is  separated  out  from  the  problem  and  universal  correction 
functions  are  available  from  the  theory  for  airfoils  at  given  incidence  and 
Mach  number  conditions. 
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FOREWORD 

This  document  constitutes  the  final  report  for  AFOSR 
Contract  F49620-80-C-0081A  for  the  period  June  1,  1980  to  May  1,  1982. 
Contained  herein  Is  a  summary  of  research  In  two  areas:  asymptotic  theory 
of  transonic  wind  tunnel  Interference  as  well  as  slender  body  theory  and 
optimization.  Each  of  the  topics  employs  an  Integrated  approach  consisting 
of  asymptotic  and  numerical  procedures  useful  In  describing  the  flow. 
Actually,  both  subjects  are  Interrelated,  since  the  three-dimensional  wind 
tunnel  problem  for  many  configurations  embodies  many  of  the  asymptotic 
Ideas  arising  In  slender  body  theory  and  the  equivalence  rule. 

This  report  Is  written  In  two  separate  sections  relevant  to  the 
previously  enumerated  areas.  These  give  a  fairly  detailed  account  of  the 
research  effort  In  each.  Principal  outcomes  of  the  program  are  the 
following: 

e  A  systematic  asymptotic  theory  for  solid  wall  Interference 
over  two-dimensional  airfoils 

e  A  new  unsteady  version  of  the  equivalence  rule 

•  A  means  of  optimizing  transonic  lifting  wing  bodies  in 
steady  flow  and  the  first  solution  our  knowledge  of  the 
axlsymmetrlc  Inverse  problem  with  closure  allowing  "nonlinear 
area  ruling" 

•  New  results  for  the  nonlinear  component  of  the  far  field 
around  lifting  airfoils  at  transonic  speeds. 

The  work  on  transonic  wind  tunnel  wall  interference  Is  a  natural  outcome 
of  activity  on  tunnel  wall  far  fields  under  previous  contracts. 

In  addition  to  the  description  of  accomplishments,  suggestions 
for  future  effort  are  provided. 
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PART  1 

SLENDER  BODY  THEORY  AND  OPTIMIZATION  PROCEDURES 
AT  TRANSONIC  SPEEDS* 

SUMMARY 

Two  areas  have  been  Investigated  within  transonic  slender  body 
theory.  In  the  first,  the  existence  of  an  unsteady  equivalence  rule  has 
been  verified  In  which  the  outer  solution  Is  governed  by  axlsymmetrlc 
unsteady  small  disturbance  theory  subject  to  an  unsteady  line  source 
Internal  boundary  condition  and  far  field.  As  In  the  steady  case,  the 
strength  of  the  line  source  depends  on  an  inner  solution  for  the  poten¬ 
tial,  which  Is  harmonic  in  cross  planes  to  the  flow.  This  solution 
specifies  the  line  source  Intensity  in  terms  of  a  rate  of  change  of  the 
cross  sectional  area  distribution  In  which  the  unsteady  motion  of  the  body 
surface  is  taken  Into  account. 

The  second  portion  of  the  research  effort  deals  with  optimization 
of  slender  lifting  wing-body  combinations.  A  restricted  class  of  these 
shapes  has  been  considered  In  which  the  wing  span  Increases  monotonlcally 
from  a  common  apex  with  the  body.  For  pointed  tails,  the  aerodynamic 
efficiency  can  be  maximized  by  minimization  of  the  wave  drag.  A  parameterized 
Inverse  (PI)  procedure  Is  utilized  to  demonstrate  substantial  decreases  In 
this  quantity.  The  PI  method  employs  a  parameterized  fairing  of  the  surface 
pressure  distribution  about  a  discontinuous  one  associated  with  a  surface 
shock.  The  parameters  are  adjusted  to  minimize  the  drag  and  satisfy  con¬ 
straints  such  as  specified  thickness,  and  base  area  as  well  as  closure  as  a 
special  case.  With  the  procedure,  the  question  of  the  existence  of  super¬ 
critical  shock-free  axlsymmetrlc  bodies  has  been  addressed.  The  examples 
Investigated  suggest  that  such  shapes  can  be  found.  Other  results  dem¬ 
onstrate  that  surface  smoothing  of  the  pressure  Is  synonomous  with  wave 
drag  elimination  due  to  focusing  and  envelope  formation  of  the  wave 
system  off  the  body. 

♦Material  from  this  section  will  be  submitted  to  the  AIAA  Journal  for 
publication  In  the  future. 
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During  the  1950's,  techniques  were  developed  to  anelyze  and  design 
aircraft  to  accelerate  through  what  was  then  called  the  sound  barrier.  The 
procedure  employed  was  an  "area  rule"1  In  which  the  drag  of  a  complete 
airplane  was  computed  from  that  of  an  equivalent  body  of  revolution. 

A  related  development  with  the  Oswatltsch  equivalence  rule2  which  was 
later  further  enunciated  by  Sprelter  and  Heaslet  In  Ref.  3.  Therein,  cases 
were  considered  In  which  the  Incidence  was  small  with  respect  to  the 
thickness.  This  approximation  leads  to  a  well  known  further  decomposition 
for  transonic  flows  which  Is  Indicated  schematically  In  Fig.  1.  As  shown, 
the  nonlinear  embodiment  Involves  the  solution  of  two  problems.  The  first 
Is  the  nonlinear  flow  over  an  equivalent  body  of  revolution  and  the  second, 
a  cross  flow  problem  which  can  In  turn  be  decomposed  Into  a  steady  asym¬ 
metric  cross  flow  and  an  unsteady  symnetrlc  flow  associated  with  the  cross 
section  expanding  at  the  rate  at  which  the  flow  moves  downstream. 

In  spite  of  Its  well  developed  status,  the  Jones -Whitcomb  area 
rule  In  which  the  equivalent  body  of  revolution  problem  Is  linear,  rather 
than  nonlinear,  can  lead  to  suboptlmal  design  solutions.  Despite  the  rela¬ 
tive  weakness  of  shocks  occurring  In  mixed  flows,  small  drag  changes  carry 
Important  mission  Implications.  An  accurate  means  of  drag  estimation  as 
well  as  a  systematic  optimization  procedure  Is  required,  particularly,  when 
the  assumptions  of  Prandtl-Glauert  theory  are  violated. 

As  one  method  of  handling  nonlinear  effects,  constrained  mini¬ 
mization  has  provided  a  useful  means  of  optimizing  transonic  shapes  under 
certain  circumstances.  In  this  procedure,  the  body  shape  Is  expressed  In 
terms  of  a  finite  number  of  parameters.  The  nonlinear  problem  for  the 
optimization  of  some  aerodynamic  quantity  such  as  drag  or  aerodynamic  effl- 
clence  as  a  functional  of  these  parameters  Is  solved  In  which  an  extrema* 
Is  sought  In  the  space  of  these  variables  using  decision  theory.  Unfor¬ 
tunately,  the  computational  effort  required  to  obtain  such  a  solution  Is 
great,  and  only  a  local  rather  titan  global  optimum  Is  usually  obtained 
In  this  manner. 
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Fig.  1.  Transonic  equivalence  rule 
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Because  of  these  deficiencies,  a  need  exists  for  other  methods. 

One  group  of  interest  deals  directly  with  desirable  loading  features  during 
the  optimization  process.  A  typical  example  occurs  in  transonic  airfoil 
theory  In  which  nearly  shockless  surface  pressure  distributions  and  mini¬ 
mization  of  adverse  pressure  gradients  are  obtained  by  a  technique  known 
as  the  inverse  method.  At  the  Science  Center,  the  approach  has  reached  a 
high  degree  of  refinement  for  airfoil  and  wing  problems.  An  account  of 
its  current  status  Is  given  In  Ref.  4.  Application  to  laminarization 
is  discussed  therein  as  well  as  the  aspect  of  closure. 

In  the  design  or  inverse  procedure,  shapes  are  sought  supporting 
preassigned  pressure  distributions  in  contrast  to  the  direct  or  analysis 
method  where  the  reverse  is  true.  A  typical  implementation  is  to  smooth 
out  jump  discontinuities  in  the  surface  pressure  distribution  that  would 
be  associated  with  shocks.  In  principle,  this  process  could  lead  to  a 
reduction  in  wave  drag.  However,  it  is  not  obvious  that  the  elimination 
of  such  a  surface  discontinuity  will  be  synonomous  with  removal  of  the 
shocks  off  the  body.  In  this  respect,  constrained  minimization  has  an 
advantage,  since  it  deals  directly  with  a  concomitant  of  the  shock,  the 
wave  drag.  We  believe  that  further  understanding  of  the  envelope  forma¬ 
tion  process  leading  to  this  discontinuity,  in  concert  with  the  inverse 
method,  could  overcome  this  deficiency  and  provide  it  with  even  greater 
power  to  solve  the  optimization  problem. 

On  the  basis  of  the  equivalence  rule  and  the  inverse  capability, 
an  attractive  possibility  occurs  regarding  optimization  of  slender  air¬ 
planes.  This  notion  has  stimulated  the  application  of  the  technology  to 
the  equivalent  body  of  revolution  component  indicated  in  Fig.  1  in  the 
research  effort  to  be  described.  Since  the  wave  drag  arises  solely  from 
this  contribution,  a  potent  means  of  improving  aerodynamic  performance 
is  available  with  simplified  versions  of  modern  computational  methods. 

Here,  the  reduction  In  complexity  arises  from  the  dominant  axial  symmetry 
of  the  asymmetric  slender  problem  In  the  nonlinear  part  of  the  flow. 
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In  what  follows,  an  account  of  the  research  concerning  the 
solution  of  the  inverse  problem  for  slender  configurations  is  given.  Prior 
to  this  discussion,  a  description  of  an  unsteady  equivalence  rule  will  also 
be  provided.  Herein, the  basic  methods  developed  in  Ref.  5  are  generalized 
to  an  unsteady  context.  Future  applications  for  the  modeling  of  flows 
relevant  to  aerodynamic  applications  are  indicated  In  this  connection.  For 
the  steady  portion  of  the  activities,  the  application  of  the  optimization 
procedures  to  achieve  minimum  wave  drag  are  described  as  well  as  salient 
Issues  dealing  with  constraints.  In  addition,  the  Implication  of  wave 
drag  reduction  measures  relative  to  aerodynamic  efficiency  is  also  given 
pointing  to  the  need  for  further  efforts  which  are  also  indicated. 
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2.0  ANALYSES 


2.1  Unsteady  Equivalence  Rule 
Basic  Equations 


Denoting  q  the  vector  velocity  p  the  density,  P,  the  pressure, 

T  the  time  and  ^  the  substantial  derivative,  the  continuity  and  momentum 
equations  for  inviscid  flow  are  respectively 


V*pq  *  0 

(la) 

3SV 

II 

II 

-1 

(lb) 

where  subscripts  denote  differentiation,  q2  =  |q)2  and  irrotationality 
has  been  assumed.  If  the  flow  is  considered  to  be  constant  state  and 
adiabatic  in  the  freestream,  the  energy  invariant  in  isentropic  flow  is 


r 


/(T) 


.2 

u2  .  a» 

"T*  y=l 


=  C2 


(lc) 


where  Is  the  velocity  potential,  U  is  the  freestream  speed,  a  Is 
the  speed  of  sound,  subscript  »  denotes  a  freestream  quantity,  and 
C  and  f  are  Independent  of  T  because  of  the  freestream  conditions 
assumed.  Substituting  (lc)  Into  (lb)  gives 

a2V-q  -  *  *n  +  ($£)  (lc*) 


where 


q2 


t 


for  cylindrical  coordinates,  in  which  R.e  are  polar  coordinates 


Rockwell  International 


SC5267.3FR 


in  planes  perpendicular  to  the  X  axis  which  Is  aligned  with  the  free- 
stream.  From  (lc‘)  the  full  potential  equation  Is  then 


(a!-*xhx  *  (,!-*r)*rr  * 


2Wex  2WeR 


(>!‘  ~fa~ ) 


'aa 

"1^  -  2  Wrx 


R2  -  2*r4>tr  -  2$x$tx  -  p-  *0$T0 


2*e*R 


-*TT=0  . 


Transonic  Unsteady  Small  Disturbance  Theory 

Heuristically  asserting  that  the  asymptotic  structure  is 
similar  to  that  of  steady  flow,  an  outer  expansion  representing  the 
flow  at  moderate  distances  from  the  slender  body  can  be  written  as 


»  >  u{x+c£a+(xtM,{)  +  ...  J 


where  c  Is  the  body  length,  6  is  the  body  thickness  ratio,  and 


1-M2 

x  *  b  ?  *  t* 1  s  K  *  “6^ fixed  as  6  ■” 0  •  (»> 


In  the  outer  limit  for  validity  for  (3a).  On  substitution  of  (3a) 

Into  (2),  the  unsteady  small  disturbance  equation 

[K-(y+1Hx3<|>xx  +  -yr  (f^)jR  +  -jr  m  2$x£  *  0  (4) 

Is  obtained.  As  In  the  steady  case,  the  preliminary  assertion  Is  made 
(pending  matching)  that  the  F  derivative  term  dominates  (4)  as  F  «►  0.  This 
leads  to  the  condition 


iir*!aw««U£3 


I 
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11m  =  S(x,t,N,n)  (5) 

f  ->  0 

where  N  and  Q  are  dimensionless  amplitude  and  frequency  parameters. 

These  quantities  arise  naturally  In  the  boundary  conditions  to  be  des¬ 
cribed.  Besides  (5),  additional  conditions  necessary  to  determine  the 
solution  of  (4)  involve  the  wake  and  far  field.  Since  our  focus  Is 
primarily  the  existence  of  an  unsteady  equivalence  rule,  the  latter  will 
be  only  briefly  and  qualitatively  Indicated  here.  For  a  subsonic  far 
field.  It  Is  reasonable  to  assume  that  an  unsteady  version  of  the  line 

J 

source  expression  given  in  Ref.  7  will  be  applicable.  The  condition 
on  the  wake  Is  determined  from  an  Inner  expansion  In  a  region  close  to 
x  axis.  The  need  for  this  representation  is  associated  with  a  non¬ 
uniformity  of  the  outer  expansion  associated  with  a  logarithmic 
behavior  given  by  an  Integral  form  of  (5),  l.e., 

<J>  =*  S(x,t)X.n?  +  g(x,t)  as  r2  0.  (5') 

Formalizing  this  Idea  and  following  Ref.  5,  the  appropriate  representation 
is 


4  =  U  {x+  62c$*(x,r*,e,t,K,n,N)  +  ...  |  ,  (6a) 

which  is  assumed  valid  In  an  inner  limit 

x.r*  *  r/6  ,  t  *  62t  fixed  as  5  0  .  (6b) 

On  substitution  of  (6b)  In  (2),  the  steady  cross  flow  Laplaclan 
equation  Is  obtained  as  In  the  steady  case  as 

♦?*r*  +  +  FT  *  0  .  (7) 
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where  the  (r*,9)  polar  coordinate  system  In  a  typical  cross  flow  plane  Is 
shown  In  Figs.  2  and  3a.  If  B(x,M)  «  0,  denotes  the  body  surface,  the 
tangency  conditions  for  (7)  is 

By  +  q  •  VB  ■  0  (8) 

If,  furthermore,  B  is  written  In  the  form 

B  *  R-6cF(x,0,T)  *  0  ,  (9) 

and  noting  that 

it  *  (l+<52$*  +  ...)  tx  +  6^Jr*  +  -6  ie  +  ...  ,  (10) 

where  lx.,lr  and  1Q  are  unit  vectors  in  the  indicated  directions,  substi¬ 
tution  of  (9)  and  (10)  Into  (8)  yields  on  retention  of  terms  to  o(6)  the 
relation 

♦e))^ _F  "  Fx(x,e,l)  +  62F^(x,0,t)  at  Fx(x,0,t)  ,  (11) 

(-?) 

where  r*  *  F  is  the  body  cross  sectional  boundary  in  the  cross  flow  plane 
shown  in  Fig.  3.  In  typical  cases,  the  function  F  can  be  written  as 

F  *  f(x,e,as1nnt)  *  f(x,8,as1nwt)  (12) 

fl  »  -jjj  •  0(1)  as  6  -*■  0 

where  if  w  is  the  dimensional  frequency,  u>  *  ~  ,  and  a  Is  an  amplitude  In 
units  of  5c.  Some  examples  follow: 
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Harmonic  Dilatation 

f  ■  aslnut  G(x,e)  (13a) 

where  aG  Is  cross  sectional  shape  at  u>t  »  j 


(13b) 


where  G  Is  the  cross  sectional  shape  In  the  pitching  or  plunging  frame, 
n  -  0  for  plunge  and  n  =  1  for  pitch,  and  Imaginary  values  of  the  radical 
In  (13b)  apply  to  e  values  outside  of  the  range  9j<  0  <  02  shown  In 
Fig.  3b.  Specialization  of  (8)  to  (13b)  Indicates  that  the  quantity  a  must 
be  0(1)  In  order  that  both  terms  of  the  left  hand  side  of  (11)  are  also 
0(1). 

For  motions  given  by  (12)  and  In  particular,  these  examples, 

Fx  +  62F^  *  Fx  +  NF^jj  COSflt  ,  (14) 

*Fx 

where  represents  differentiation  with  respect  to  the  third  argument, 
and  N  »  a Q62  *  o(<52)  since  a,  fi  have  been  assumed  to  be  0(1).  If  Si  ■  o(l) 
and  a  *  o(&~2),  unphysical ly  large  amplitudes  would  occur,  or  the  radial 
variation  would  dominate  the  circumferential  ones  In  the  left  side  of  (11), 
departing  from  desired  generality.  If  n=o(6-2),  with  a  *  0(1),  then  addi¬ 
tional  temporal  derivative  terms  would  be  required  in  the  equations  of 
motion  and  the  boundary  conditions  as  In  the  linear  case.  Note  that  for 
planar  flows,  (of  Eq.  (4-7)  of  Ref.  6),  the  expression  analogous  to  (14) 
doee  contain  the  second  term  involving  an  amplitude  parameter  such  a  N. 

There,  the  Inclusion  does  not  lead  to  large  amplitudes. 

Using  Identical  procedures  to  that  given  In  Ref.  5,  the  boundary  con¬ 
dition  (14)  can  be  written  as  the  following  Neumann  condition  In  the  cross  plane: 


3<t* 

3n 


1-12 


9 


Rockwell  international 

Science  Cantor 

SC5267.3FR 


This  is  explicitly  the  same  boundary  condition  (except  with  t  dependencies) 
as  for  the  steady-state  inner  solution  of  Ref.  5.  Using  Green's  theorem 
as  in  Ref.  5,  the  asymptotic  behavior  of  for  large  r*  is 


♦*  *  jz  Unr*)Ax 


(16a) 


where 


i  r2v 

•U  F!de 


(16b)* 


In  accord  with  the  previous  remarks,  the  complete  inner  solution  can  be 
written  as 

<j>*  «  S*(x,t)  log  r*  +  g*(x,t) 


♦2 

n*l 


9n(x»t)cosn0  +  fn(x,t)sinne 


„*n 


(17) 


The  matching  procedures  developed  in  Ref.  5  applfed  here  give  the  following 
results 


S(x,t)  »  S*(x,t)  ■  jn  Ax 


(18a) 


♦The  limits  shown  are  applicable  If  r*-F  encloses  the  origin  as  a  star 
shaped  region,  l.e.,  r*«F(0,x)  Is  a  single  valued  function  of  8.  If  not, 

then  A  ■  £  F2de,  where  the  0  limits  are  selected  so  that  the  contour  r*»F 

Is  traversed  In  the  counterclockwise  sense. 
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g(x,t)  =  g*(x,t)  ,  (18b) 

with  the  Inner  and  outer  expansions  given  respectively  by 

*INNER  "  U  jx+c[(62iog  6)2S(x,t)+62d»(x,r*,t;K,n)  +  ...]|  (19a) 

for  (6b)  and, 

*  OUTER  =  U  |x+c62<{i(x,r,t;K,n)  +  ...  j  (19b) 

for  (3b),  and  the  52fcog  6  switchback  term  Inserted  for  matching  purposes. 

From  the  energy  Invariant  and  Isentropy,  the  surface  pressures 
on  the  body  r*=F  are  given  by 


2Sx(x,t)  *og62F  +  2gx(x,t)  +  vj  +  Wg 


(21a) 


where 


5  “  gncosn9  +  fns1nn0 

’TLt - ^*1) - 

n-1  h 


"B'S 

n»l 


n  |  -gns1nn9  +  fncosn6 } 


In  which  g„  and  f_  are  functions  of  x  and  t. 
n  n 


(21b) 


(21c) 
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Summary  -  Unsteady  Equivalence  Rule 

The  unsteady  equivalence  rule  Is  embodied  In  Eqs.  (21)  In 

which 


9 


lim  -  S(x,t)  tog  ? 

r  0  < 


(22) 


with  S  determined  by  (18a)  and  $  the  outer  perturbation  potential  Is 
the  solution  of  the  unsteady  small  disturbance  equation  (4),  In  which 
*  0,  due  to  the  equivalent  body  of  revolution  boundary  condition  (5). 
The  wake  condition  by  symmetry  Is 

^(x.0,1)  «  0  ,  x  >  1  (23) 

In  the  solution  of  the  outer  problem. 

In  a  practical  application,  the  procedure  described  pre¬ 
viously  could  be  utilized  to  compute  transonic  unsteady  flows  over  low 
aspect  ratio  configurations  such  as  fighter  arrangements.  To  obtain 
a  truly  useful  tool  however,  generalizations  are  required  to  handle 
shapes  In  which  the  wings  are  non-monotonically  increasing  In  span 
from  the  nose  of  the  wing-body  arrangement.  For  these  cases,  the 
wakes  play  an  important  role,  modifying  the  autonomy  df  downstream 
sections.  As  in  steady  Interactions,  the  higher  order  matching  is 
Important,  including  a  theory  of  the  materialization  of  shocks  in  the 
near  field. 
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2.2  Steady  Slender  Body  Theory  and  Optimization  Studies 

The  previous  results  when  specialized  to  steady  flow  can  be  used 
to  study  L/D  optimization  of  wing-body  configurations  of  the  type  shown  In 
Fig.  4.  Here,  the  body  has  a  circular  cross  section,  and  the  wing  Is 
Increasing  monotonlcally  from  the  pointed  apex  to  the  base  section  and  Is  of 
zero  thickness.  Based  on  Ref.  5  and  8, the  aerodynamic  efficiency  at  an  angle 
of  attack  a  Is  given  as 


L  2”ts2-ft02+V/s2)a 

^  6k  Jb+I 


(24a) 


where 


B  -  trS(l)g(l) 


1 

I 


J2V(l,F,0)FFx(l,0)de  -  [log  fi][2irS2(l)]  (24b) 


I  *  -2tt 


'(x)dx 


(24c) 


Rq  «  body  base  radius  In  units  of  the  chord 


s  *  span  at  base  In  units  of  the  chord. 


and  g  and  S  are  steady  specializations  of  that  given  In  (22)  and  (18a). 

In  general,  the  specification  of  a  fixed  base  could  be  driven  by  an  engine 
packaging  requirement  or  RCS  considerations. 
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configuration  for  L/D  optimization  studies 


6, 
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In  the  research  conducted  during  the  course  of  the  program  a 
major  emphasis  was  placed  on  drag  reduction  associated  with  the  minimization 
of  the  denominator  In  (24a).  For  a  fixed  Rq/s,  (24a)  indicates  that  the 
maximum  L/D  corresponds  to  a  minimum  wave  drag.  The  minimization  Involves 
global  dependencies  that  modify  B  through  g(l)  and  4>*(1,F,6)  In  (24b).  To 
simplify  the  considerations,  a  part  of  the  Investigation  was  limited  to  a 
simpler  situation  In  which  Rq*0,  l.e*,  the  body  Is  pointed  at  Its  tall. 

For  this  case,  (24a)  specializes  to 


l.e.,  the  lift  Is  given  by  a  wing  alone  value  and  the  boundary  terms  are 
absent  from  the  drag.  For  a  fixed  span  at  the  base,  the  aerodynamic  effi¬ 
ciency  Is  again  maximum  when  the  denominator  Is  minimized.  For  (24a'), 
the  complicated  global  dependencies  In  the  boundary  terms  are  not  present. 

To  meet  the  foregoing  minimization  requirement,  an  effort  In 
developing  and  applying  special  methods  to  derive  pointed  bodies  supporting 
shockless  surface  pressure  distributions  will  be  described.  These  proced¬ 
ures  are  different  than  another  one  that  evolved  in  an  early  phase  of 
this  effort  and  was  reported  In  Ref.  9.  The  latter  had  the  focus  of  deal¬ 
ing  with  certain  coupling  In  the  boundary  conditions  peculiar  to  the  present 
axisymmetric  generalization  of  the  two-dimensional  airfoil  problem  described 
In  Ref.  4.  In  what  follows,  a  description  will  be  given  of  a  new  method 
tha1;  provides  a  treatment  of  constraints  that  are  vital  for  the  L/D  opti¬ 
mization  problem  related  to  equations  such  as  (24a)  and  (24a‘),  and  In 
addition  gives  the  first  method  of  obtaining  closure  to  our  knowledge  for 
the  axisymmetric  inverse  problem. 
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This  "parametrized  Inverse"  (PI)  procedure  consists  of  the 
following  steps; 

1.  Starting  with  a  known  body*  the  analysis  boundary  value 
problem  consisting  of  Eqs.  (4),  (5),  (18a)  and  (23)  is  solved 
to  obtain  a  surface  pressure  Cp  given  by  (21). 

2.  In  general,  the  pressure  distribution  of  Step  1  will  have  a 
discontinuity  associated  with  a  shock  terminating  the  super¬ 
sonic  zone.  A  fourth  order  polynomial  is  used  to  "refair" 
the  region  in  the  vicinity  in  the  shock  as  schematically 
Indicated  In  Fig.  5.  In  the  notation  of  the  figure,  control 
points  at  Xj  and  X2  are  shown,  and  continuity  of  slope  and 
value  of  Cp  Is  required  at  these  locations.  Two  parameters 
are  Introduced  toward  minimizing  the  drag  Cq  and  obtaining 
closure  in  the  sense  that  F(1)*0.  (A  more  general  case  In 
which  F(l)=Fj  is  prescribed  is  also  discussed  in  what  follows.) 
The  first  parameter  a  is  the  midpoint  value  of  Cp  shown  In 
Fig.  5,  l.e. ,  Cp(x3).  The  second,  is  the  intensity  y  of  the 
leading  edge  singularity,  where  for  bodies  pointed  at  tha 
nose  and  tall,  according  to  (21a), 

CD  *  yfcog(l±x)  as  x-*-+l,  (25) 

Kb  (If  F~l±x) 

and  the  subscript  b  refers  to  the  body  value.  A  Cp  variation 
near  the  nose  of  a  symmetric  parabolic  arc  body  illustrating 
this  behavior  Is  shown  In  Fig.  6,  where  the  freestream  Mach  t 
number  M*  *  0.98  and  the  thickness  ratio  6  ■  0.167.  It  can  be 
seen  that  the  numerical  procedure  tracks  the  singularity  quite 
well,  except  the  first  mesh  point,  which  Is  understandable. 

In  fact  the  proper  strength  y  Is  also  obtained  on  Inspection  of 
*  the  slope  of  the  curves  and  (5*). 


9.  * 
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r 


c|>  <x3> E* 


Fourth  order  smoothing  used  for  refairing 
shocked  surface  pressure  distributions 
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1+x 


I 

Fig.  6. 


Singular  behavior  at  the  nose  of  a  fore  and  aft 
symnetrlc  parabolic  body  H»  ■  0.98,  6  ■  0.167 
(from  SLOR  solution) 


Rockwell  International 

leiemce  Cantor 


SC5267.3FR 


3.  In  the  Iteration  procedure,  a  given  value  of  the  parameter  a 
Is  assumed  as  well  as  y.  The  latter  Is  regarded  as  a  constant 
multiplicative  factor  for  the  entire  Cp  (x,a,y)  distribution 
on  the  Interval  -l<x<l.  b 

Regarding  (21a)  as  an  ordinary  differential  equation  for  F  in 
which  gx  Is  assumed  known  from  a  previous  Iteration,  with  Cp 
specified,  Wg»0,  and  Vg«F',  a  two  point  boundary  value  problem 
Is  solved  for  this  equation  In  which 

F(0)  *  F(l)  «  0  (26) 


are  the  boundary  conditions.  With  gx  fixed  at  this  stage, 
the  solution  technique  Is  the  method  of  bisection  coupled 
with  a  Runge-Kutta  Integration  procedure.  In  the  Implemen¬ 
tation,  two  Initial  values  of  y,  Vj,  and  v>2  are  obtained  to 
bound  a  range  where  R(l)  changes  sign.  Values  of  y  are  changed 
in  the  left  hand  side  of  (21a)  using  relaxation.  A  set  of 
loops  are  Initiated  In  which 


cp  «=  rcp  +  (l-r)  Cp 

"RELAXED  r0LD  "TARGET 


(27a) 


and 


C 


P 


NEWER 


(27b) 


The  ordinary  differential  equation  (21a)  Is  solved  with 
(27b)  and  the  bisection  method  until 


|F(1) |<.001. 


(27c) 
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Our  experience  has  been  that  convergence  Is  very  rapid  with 
this  approach.  Typical  values  for  the  yj  are:  Uj  *  0.2  and 

yg  *  1.8. 

4.  Once  a  new  F  distribution  Is  obtained,  the  analysis  algorithm 
described  In  Step  1  Is  used  to  recalculate  Cp  and  g'(x).  In 
this  step,  F  is  relaxed  In  accord  with  the  relationship. 


_  fold  +  fnew 

RELAXED  2 - 


(28) 


5.  Steps  3-5  are  repeated  until 

Cp  *  Cp 

*NEW  ^TARGET 

6.  The  value  of  the  drag  coefficient  Is  correlated  against  other 
values  of  the  parameter  a  Introduced  In  Step  3.  In  (27a),  the 
value  of  r  Is  adjusted  In  accord  with  the  rule 

r  *  !-  ,  (29) 

where,  l<n<40,  and  n  *  number  of  Iterations.  Thus,  according 
(29),  T-*0  as  n-*«. 

7.  In  a  recent  development,  the  value  of  F  Is  renormalized  by 
Introducing  a  new  F  In  accord  with  the  relation 

F  «  yJ-  (30) 

rm 

Step  (7)  Is  Integrated  Into  the  Iteration  procedure  as  a  means 
of  Introducing  a  fair  basis  of  comparison  of  the  drag  levels. 

It  should  be  noted  that  In  the  SLOR  analysis  code  employed  In  Step  4, 

It  was  adequate  to  use  only  15  sweeps  before  the  next  phase  of  the  Iteration 
process  was  Initiated. 
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In  what  follows,  Illustrations  of  the  PI  method  previously  described 
to  obtain  shockless  pressure  distributions  and  closed  axl symmetric  bodies 
will  be  provided.  The  Implications  of  these  measures  with  respect  to  wave 
drag  will  also  be  addressed.  Also  to  be  Indicated  are  calculations  from 
algorithms  devised  for  the  satisfaction  of  other  constraints  besides  closure. 
Viscous  effects  and  base  drag  on  finite  base  afterbodies  have  not  been 
addressed  here  but  deserve  consideration  In  future  effort.  Also,  basic 
Issues  of  existence  and  uniqueness  deserve  attention.  Another  aspect  Inves¬ 
tigated  was  the  attainability  of  a  zero  wave  drag  body  of  revolution. 

To  show  the  capability  of  the  Inverse  method  to  eliminate  surface  shocks, 
modifications  of  symmetric  parabolic  arc  bodies  were  considered  Initially.  Fig.  7 
depicts  a  vertically  expanded  version  of  such  a  body  (dashed  curve)  where  the 
notation  of  R  Is  used  In  the  figures  for  F.  The  associated  pressure  (dashed 
curve)  Indicates  a  jump  discontinuity  produced  by  a  shock  at  a  location  downstream 
of  the  midpolntof  the  body.  For  these  calculations,  89  points  were  used  In  the 
x  direction  with  55  points  across  the  body.  In  the  r  direction,  31  points  were 
employed.  Clustering  was  used  in  the  vicinity  of  the  nose  and  tall  of  the  body 
for  the  x  grid  and  the  usual  expansion  of  the  r  grid  at  large  distances  from  the 
axis  of  symmetry  was  also  applied.  A  smoothed  pressure  distribution  shown  as 
the  solid  curve  arises  from  the  4th  order  polynomial  discussed  previously,  where 
the  hump  Is  indicative  of  the  solution  for  the  parameter  a  required  to  achieve 
closure.  The  closure  of  the  redesigned  body  shown  as  the  solid  curve  Is  typical 
of  the  success  we  have  had  with  the  PI  procedure  In  meeting  this  goal. 

In  Fig.  8,  a  schematic  of  the  formation  of  a  shock  over  an  airfoil 
Is  shown  In  which  the  discontinuity  Is  created  by  reflected  waves  of  second 
family  characteristics  from  the  sonic  line. 

Fig.  9  Indicates  a  typical  wave  system  over  a  body  of  revolution. 

In  connection  with  these  figures,  It  Is  of  Interest  to  establish  under  what 
conditions  envelopes  form,  and  more  Importantly,  when  removal  of  the  surface 
pressure  discontinuity  results  In  zero  wave  drag.  It  should  be  noted  In  this 
connection,  that  solutions  of  Ref.  9  as  well  as  the  PI  method  of  this  report 
produce  reduced,  but  finite  wave  drag  with  smoothed  surface  pressure  distri¬ 
butions.  In  Ref.  9,  the  Issue  of  closure  was  not  addressed.  Subcrltlcal 
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Fig.  7.  Body  redesign:  H  *  0.98,  6  ■  0.167 
K  -  (1-Hl)/M&62  ■  1.48,  (Run  40H) 
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Fig.  8.  Envelope  formation  of  a  shock  over  an  airfoil 


Fig.  9.  Shock  structure  over  a  hotly  of  revolution 
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cases  discussed  therein  demonstrated  a  finite  drag  for  unclosed  bodies. 
Since  a  finite  base  body  violates  the  assumptions  In  the  proof  of 
D'Alembert's  paradox,  this  result  should  not  be  surprising.  Because  the 
bodies  can  be  closed  In  the  PI  method,  the  wave  drag  can  be  traced  directly 
to  the  shock  system. 


Regarding  the  Issue  of  drag,  the  basic  drag  coefficient  Cp  of  the 
parabolic  body  for  the  case  of  Fig.  7  was  computed  to  be  .0459.  For  the 
refaired  pressure  distribution  obtained  In  the  Iterations  of  the  PI  method, 
convergence  to  a  value  of  a— 0.275  lead  to  Cp  ■  0.00022,  practically 
zero  to  within  the  truncation  error  of  the  calculations.  This  result  Is 
tempered  by  the  reduction  of  R^  from  unity  for  the  parabolic  case,  to  a 
value  of  0.524  In  the  refaired  case.  Since  there  Is  a  natural  reduction  In 
drag  with  thickness  ratio,  a  fair  comparison  between  the  original  and 
designed  body  should  really  be  performed  with  both  at  the  same  thickness. 

The  drag  of  the  rescaled  redesign  Is  obtained  using  the  transonic  similarity 
rule,  which  for  axl symmetric  bodies  of  thickness  ratio  6  Is 


_  CD 


CD  S  6*  *  f<K 


Ola) 


where 


K'»  (l-Mw) /62. 


(31b) 


Accordingly,  If  a  CD  Is  associated  with  an  Ri  =  p,  and  p<l,  say,  then 

MAX 

the  body  thickness  ratio  should  be  Increased  by  a  factor  p"1 


coefficient  C, 


9? 


Is  associated  with  an  Ri  =  p, 

and  a  new  drag 

should  be  calculated  from  a  new  analysis  solution  to  evaluate 


f(IC,),  where  Kg  *  p2K^.  Then 


(32) 


is  the  appropriate  comparison  drag  coefficient  for  the  blown-up  redesigned 
body*  Since  the  rescaling  of  Kj  to  p2Kj  drastically  changes  the  flow  field. 
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with  a  probable  generation  or  Intensification  of  shocks  for  y<l,  It  Is 
Important  to  Implement  Step  7  In  the  PI  procedure  previously  described  If 
a  meaningful  drag  reduction  Is  to  be  achieved.  Some  results  of  preliminary 
efforts  In  this  direction  will  be  described  subsequently.  The  procedure 
described  In  Ref.  9  embodies  the  feature  of  renormalizing  R^  during  Its 
Iteration  cycle. 

As  a  basis  of  a  different  sort  of  standardized  comparison,  the 
flow  field  over  a  reference  parabolic  arc  body  with  R^  ■  .524  was  com¬ 
puted  at  K  *  1.48,  the  value  utilized  for  the  original  parabolic  and  the 
redesigned  body.  By  contrast  to  the  redesigned  body,  a  surface  shock 
was  exhibited  for  the  parabolic  body,  although  the  drag  was  still  low 
(CD  *  .00066),  but  three  times  larger  than  the  redesigned  body. 

Since  the  emphasis  In  the  PI  method  Is  to  achieve  closure  as  well 
as  address  the  question  of  the  existence  of  drag-free  bodies  at  transonic 
speeds,  the  Issue  of  standardization  Is  of  lesser  Importance  at  this 
stage  of  the  discussion.  The  significance  of  the  redesign  of  Fig.  7  can  be 
more  fully  appreciated  from  a  picture  of  the  flow  above  the  body.  These  are 
shown  In  what  follows  In  terms  of  Isobar  and  IsoMach  level  lines.  In  Fig.  10, 
the  Isobar  pattern  for  the  parabolic  reference  body  at  the  conditions  of 
Fig.  7  are  Indicated.  The  solid  lines  indicate  expansion  regions  and  the 
dashed  lines  signify  compressions.  Evident  In  the  map  is  the  concentration 
of  contours  above  the  75%  location  of  the  body  associated  with  the  shock. 

A  corresponding  IsoMach  map  Is  shown  In  Fig.  11,  where  the  supersonic  bubble 
Is  shown,  as  well  as  the  transition  from  supersonic  to  subsonic  flow  at 
the  shock. 

By  contrast  to  the  shocked  pattern  evident  In  Figs.  10  and  11, 

Figs.  12  and  13  Indicate  a  shock-free  mixed  flow.  Thus,  In  spite  of  the 
Rpi^  reduction,  these  figures  provide  some  Insight  Into  the  Issue  raised 
previously  regarding  whether  drag-free  bodies  of  revolution  exist  in  mixed 
flow.  On  the  basis  of  the  apparent  shock-free  nature  of  the  flow  field 
and  the  extremely  low  drag  value  CQ  ■  .00022,  prospects  would  appear 
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Fig.  12.  Isobar  map  for  redesigned  body 
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Fig.  13.  IsoMachs  -  redesigned  body,  (height  of  sonic  bubble  *  .24) 
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encouraging  for  finding  such  bodies.  The  PI  method  appears  to  provide 
an  attractive  procedure  to  generate  such  shapes  that  close.  To 
further  elucidate  the  drag  controlling  mechanisms  relevant  to  this 
evidence,  the  nonlinear  axlsymmetrlc  hodograph  should  receive  attention. 

In  two-dimensional  airfoil  flows,  the  hodograph  is  linear,  Garabedlan 
solves  It  numerically,  and  Nleuwland  provides  analytical  solutions. 

In  some  sense,  these  treatments  prove  the  existence  of  shock-free  airfoils. 
Due  to  the  nonlinearity  of  the  hodograph  for  the  axlsymmetrlc  case, 
existence  of  shock-free  bodies  can  probably  only  be  demonstrated  In  a 
constructive  sense,  l.e.,  the  numerical  solution  of  the  boundary  value 
problem  In  the  hodograph  plane. 

It  Is  also  interesting  to  note  that  Fig.  7,  shows  that  the  point 
of  maximum  thickness  moves  rearward  for  drag  reduction  In  mixed  flow 
over  a  body.  This  Is  in  agreement  with  the  optimum  supersonic  airfoil, 
which  Is  a  vertically  symmetric  wedge  section,  whose  maximum  thickness 
moves  aft  with  increasing  Mach  number. 

Returning  to  Fig.  5,  a  study  of  the  sensitivity  of  the  drag 
due  to  the  fairing  parameter  a  was  conducted.  The  results  are  shown 
in  Fig.  14,  where  CD  and  are  plotted  versus  -a.  A  systematic  trend 
is  Indicated  showing  that  potent  gains  can  be  achieved  through  the  use 
of  relatively  simple  minded  smoothings  of  the  shock  discontinuities. 

In  a  previous  treatment  of  the  optimization  problem  for  bodies 
of  revolution,  Chan19  treated  a  two-parameter  family  of  general  parabolic 
bodies  of  the  form  S(x)  *  Axn(l-x)m,  A  ■  (n*m)^n4mV8nlWn  using  "numerical 
optimization"  to  determine  the  best  selection  of  n  and  m  to  minimize  the 
drag.  As  In  most  optimization  procedures,  the  best  selection  of  n  and 
m  did  not  eliminate  the  drag,  and  provided  only  a  local  extremum.  In  the 
PI  procedure,  there  appears  to  be  an  opportunity  to  control  the  drag 
directly,  through  the  shock  smoothing  mechanism  mentioned  previously. 


1-34 


Rockwell  International 


PCffltf  wtHHf 

SC5267.3FR 


Fig.  14.  Nave  drag  and  EBR  maximum  thickness  Ry.x  versus 
surface  CD  smoothing  parameter,  a  ** 
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Another  study  of  wave  drag  associated  with  various  values  of  the 
parameter  a  Is  shown  In  Fig.  15,  where  the  entropy  jump  integral  across 
the  shock  which  Is  proportional  to  the  wave  drag  Is  shown.  Dramatic  differ¬ 
ences  are  shown.  In  which  a  fairing  parameter  a  *  -.109  can  be  selected  to 
give  a  drag  actually  higher  than  that  of  the  parabolic  body.  Also  Indicated 
is  a  very  small  value  associated  with  a  *  -.2.  In  the  figure,  the  Irregu¬ 
larities  Indicated  are  associated  with  the  usual  difficulties  of  shock  capture. 

In  the  study  of  the  optimization  process,  a  question  of  significance 
arises  concerning  whether  the  parabolic  body  is  a  disadvantageous  Initial 
iterate  for  PI  process.  The  thrust  of  the  question  pertains  to  the  pos.  1- 
bllity  that  flattened  bodies  such  as  supercritical  airfoils  are  useful  as 
drag  reducers,  since  their  slope  contributions  to  the  wave  drag  Integral  are 
nearly  zero  on  their  nearly  cylindrical  portions.  The  associated  shocks 
could  also  be  weak  due  to  suppression  of  the  envelope  of  the  reflections.  One 
candidate  for  consideration  Is  the  shape  shown  In  Fig.  16  hereinafter  referred 
to  as  the  tanh  body.  In  spite  of  the  deleterious  effect  of  overexpansions 
near  the  forebody  shoulder  and  the  resulting  shock,  the  drag  can  be  mini¬ 
mized  to  about  1/2  of  the  parabolic  value  through  the  selection  of  the 
parameter  e,  the  effective  forebody  length.  It  Is  interesting  to  note  that 
e  *  0.12  gives  about  three  times  the  parabolic  body  result  for  the  tanh  shape. 

As  an  Illustration  of  some  of  the  nonlinear  wave  steepening  mech¬ 
anisms  present.  Fig.  17  shows  a  sequence  of  axial  distributions  of  pressure 
at  various  heights  above  the  body.  The  results  are  plotted  for  the  minimum 
drag  value  of  e.  The  Impossibility  of  Improvement  stems  from  overtaking  of 
the  expansion-compression  hump  over  the  afterbody  shoulder  by  another 
hump  from  the  forebody  shoulder.  As  the  height  Increases,  the  latter 
swallows  up  the  former.  At  even  greater  heights,  a  mild  shock  forms. 

The  shock  was  detected  In  this  case  by  monitoring  finite  jump  discontin¬ 
uities  of  $xx  and  using  expanded  scales  In  graphic  visualizations  of  the 
process. 
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Fig.  15.  Entropy  jump  shock  drag  for  various  bodies 
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Fig.  17.  Development  of  Cp  vs.  x  at  various  levels  above  tanh  body 
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To  demonstrate  the  Implementation  of  Step  7  In  the  PI  method,  the 
case  of  Fig.  7  was  run  a-galn  In  which  the  parameter  a  ■  -0.17.  In  Fig.  18, 
the  body  solution  Is  shown  as  the  Indicated  bulge  about  the  parabolic  con¬ 
tour.  Fig.  19  Indicates  the  redesigned  pressure  distribution  associated 
with  the  bulged  contour.  It  should  be  noted  that  spikes  are  due  to  limited 
convergence  of  the  iterations  for  this  case.  Part  of  the  marginal Ity  of 
convergence  stems  from  possible  existence  and  nonexistence  of  neighboring 
solutions  In  the  Iteration  space.  However,  these  difficulties  are  deemed 
to  be  solvable  with  further  numerical  experimentation.  The  axial  distri¬ 
butions  of  pressure  at  various  heights  above  the  body  have  the  same  steep¬ 
ening  characteristics  that  are  Indicated  In  Fig.  17. 

The  Important  point  about  the  cases  In  Figs.  17  and  18  Is  that 
they  conclusively  demonstrate  that  smoothing  on  the  body  may  not  always 
eliminate  wave  drag  or  shocks  off  the  body.  Therefore,  the  envelope 
formation  process  should  be  studied  carefully  in  future  effort. 

In  connection  with  earlier  remarks  on  L/D,  the  last  example  shown 
In  Fig.  20  Illustrates  the  use  of  the  PI  procedure  for  a  constrained  finite 
base.  As  usual,  the  reference  parabola  Is  indicated  for  comparison  purposes. 

The  associated  refaired  pressure  distribution  is  Indicated  In  Fig.  21.  For 
this  case,  the  maximum  velocity  potential  differences  between  old  and  new 
solutions  near  the  end  of  the  Iteration  sequence  employed  were  as  follows: 

{  ••  Nax  (♦"+1  -  ♦")  -  0.0003 
:  Max  ($n+1  -  4,n)  -  0.00815 


Design 

Loop 


Fig.  19.  Pressure  distribution  (partially  converged)  for  redesigned 
body  with  R|^  ■  1,  case  of  Fig.  18 
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Body  design  with  a  constrained  finite  base 
parameters  Identical  to  Fig.  18 
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The  convergence  rate  of  the  second  difference  was  slow,  but  the  Impact  of 
this  behavior  appears  to  be  small  on  the  difference  between  the  actual  value  of 
Cp  obtained  and  the  target  variation.  In  Table  1,  a  comparison  of  the  con¬ 
strained  and  RgASE  cases  is  shown.  Drag  levels  are  also  indicated 
against  the  parabola  value. 

TABLE  1.  COMPARISON  OF  DRAGS  -  OPTIMIZED  CONSTRAINED  BODIES 
AGAINST  PARABOLIC  TYPE 


Run 

a 

rbase 

^X 

CD 

cd/rmax 

90A 

-0.17 

0.2 

0.841 

0.0130 

0.026 

90B 

-0.17 

0 

1 

0.0384 

0.0384 

PARABOLA 

- 

0 

1 

0.0459 

0.0459 

Jt  is  evident  from  the  table  that  substantial  improvement  from  the  drag 
level  of  the  parabolic  body  was  made  with  the  constrained  maximum  thickness 
generalization  of  the  PI  procedure.  In  the  implementation  for  the  finite 
base  case,  the  absence  of  the  tail  singularity  was  Ignored,  with  no  major 
effect  on  the  convergence  of  the  closure  Iteration  loop.  Worthy  of  mention 
in  this  connection  is  the  subtle  distinction  between  a  sting  support,  open 
body,  wake  and  finite  base.  For  the  first,  second  and  fourth  situations, 
a  tail  singularity  is  possible,  whereas  for  the  third,  it  is  unlikely. 
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4.0  CONCLUSIONS 

In  the  previous  sections,  results  of  oar  research  on  transonic 
slender  body  theory  toward  obtaining  an  unsteady  equivalence  rule  and 
application  of  an  axisynmetric  Inverse  design  method  to  the  study  of 
nonlinear  area  ruling  concepts  have  been  presented.  Asymmetric  configura¬ 
tions  were  considered  In  the  research.  For  certain  classes  of  these 
shapes  In  which  the  wing  Is  of  zero  thickness  and  Increases  its  span 
montonically  from  a  cannon  apex,  the  L/D  optimization  problem  consists 
of  minimizing  the  drag  of  the  body  for  a  fixed  span  of  the  wing. 
Furthermore,  if  the  body  is  pointed,  the  wing-body  Interference  lift  has 
been  shown  by  previous  Investigators  to  be  that  of  the  isolated  wing. 

For  the  drag  minimization  process  for  these  lifting  wing  body  combinations, 
it  is  clear  now  that  the  PI  method  developed  during  the  research  program 
can  provide  an  extremely  useful  tool. 

Substantial  reductions  of  the  wave  drag  can  be  achieved  In  many 
cases  In  which  the  surface  pressure  jump  discontinuity  is  refaired.  It 
has  been  also  demonstrated  that  the  refairing  process  can  be  utilized  to 
achieve  satisfaction  of  constraints  such  as  closure  and  fixed  maximum 
thickness.  However,  there  are  certain  situations  in  which  shocks  develop 
off  the  body.  These  have  been  demonstrated  In  our  effort.  For  these 
cases,  the  elimination  of  a  surface  discontinuity  may  reduce  the  wave  drag 
but  not  completely  eliminate  It.  Therefore,  a  knowledge  of  the  theory 
of  characteristic  propagation  and  envelope  formation  for  flow  over  axl- 
symmetrlc  bodies  could  provide  target  surface  pressures  to  eliminate 
these  off-body  shocks.  Other  aspects  of  the  research  have  shown  that 
the  Initial  Iterate  In  the  design  process  can  play  a  role  In  the  degree 
of  drag  minimization  that  is  possible.  In  fact,  a  study  of  a  nearly 
Cylindrical  body  when  optimized  In  Its  forebody  and  afterbody  sections 
provides  substantial  reductions  from  the  symmetric  parabolic  shape,  if 
the  lengths  of  these  sections  are  optimized.  The  role  of  overexpansions 
and  shock  formation  Is  Intimately  Involved  In  the  proper  adjustment 
of  these  lengths  as  has  been  previously  discussed. 
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One  question  answered  by  the  research  effort  concerns  the 
existence  of  shock-free  flow  over  axlsymmetrlc  shapes.  From  the  results. 
It  would  appear  that  to  within  the  truncation  errors  of  the  numerical 
methods,  the  answer  must  be  in  the  affirmative.  Because  of  the  limited 
Information  available  however,  more  cases  should  be  Investigated  and  the 
existence  question  should  be  studied  by  other  methods  such  as  the  axl¬ 
symmetrlc  hodograph. 

For  unsteady  flows,  analytical  results  have  been  provided  that 
Indicate  the  existence  of  an  unsteady  equivalence  rule.  Oust  as  In  the 
steady  case,  the  three-dimensional  nonlinear  problem  for  slender  shapes 
can  be  simplified  to  an  axlsymnetrlc  one,  with  Internal  boundary  condi¬ 
tions  associated  with  asymptotic  behavior  of  an  Inner  solution  which 
Is  harmonic  In  cross  planes.  With  the  substantial  current  cost  of 
computing  transonic  unsteady  flows,  as  well  as  the  need  for  avoiding 
non-classical  flutter,  the  simplifications  and  cost  reductions  afforded 
by  such  a  rule  could  be  of  substantial  advantage. 
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5.0  RECOMMENDATIONS 

On  the  basis  of  the  previous  findings.  It  Is  evident  that 
slender  body  theory  may  be  exploited  with  modern  computational  techniques 
to  provide  a  basis  for  design  synthesis  of  three-dimensional  fighter 
and  missile  configurations.  Current  practice  In  many  organizations  Is 
to  use  a  linear  "base  point"  and  optimize  on  a  trial  and  error  basis 
from  direct  solvers.  Inverse  methods  applied  to  an  Isolated  component 
such  as  a  wing,  or  wind  tunnel  tests.  Because  this  may  result  In  a 
suboptlmal  extremum  prejudiced  by  the  linear  base  point,  an  alternate 
procedure  Is  needed.  For  this  purpose,  the  PI  method  described  in  this 
report  appears  attractive.  However,  for  this  tool  to  be  more  useful, 
additional  development  is  required  along  the  following  lines: 

1.  Wings  of  zero  thickness  In  combination  with  circular 
arc-cross  section  bodies  were  investigated.  Although  pro¬ 
viding  an  indication  of  the  significance  of  wing  body 
Interactions  as  they  relate  to  appropriate  directions  for 
gains,  the  Issues  of  wing  thickness  and  body  asymmetries 
should  be  addressed.  For  such  generalizations,  an  Involved 
tradeoff  becomes  evident  between  exposed  area,  for  lift 
and  nonlinear  area  ruling  required  to  reduce  the  wave 
drag.  An  associated  challenge  Is  related  to  the  ability 

to  "unfold"  the  optimum  equivalent  body  of  revolution 
determined  by  the  PI  method  to  produce  optimum  compromises 
of  design  attributes  such  as  those  between  aerodynamic 
efficiency  and  volume. 

2.  Another  aspect  of  significance  In  the  optimization  process 
deals  with  the  Interaction  of  wakes  associated  with  plan- 
forms  of  non-montonlcally  Increasing  span  In  the  downstream 
direction.  The  Implication  of  these  effects  on  the  equiva¬ 
lence  rule  and  the  optimization  procedures  considered  here 
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should  be  considered  In  future  work.  Propulsive  effects 
such  as  jets  emanating  from  the  base  are  other  possible 
avenues  of  research  In  which  streamwlse  Interactions  are 
of  Importance.  Effort  should  be  expended  on  the  exploita¬ 
tion  of  these  effects  for  beneficial  Interference. 

3.  On  a  more  fundamental  basis*  the  theory  of  characteristic 
propagation  should  be  utilized  to  shed  light  on  envelope, 
focusing  and  caustic  phenomena  associated  with  shock  wave 
development  off  axlsymmetrlc  bodies,  particularly  when  no 
surface  discontinuities  are  encountered  as  In  the  examples 
described  in  previous  sections.  This  Information  could  be 
used  to  suggest  the  types  of  faired  pressure  distributions 
leading  to  zero  wave  drag.  Although  we  have  demonstrated 
encouraging  evidence  regarding  the  existence  of  shock-free 
axi symmetric  shapes,  more  information  is  needed  on  how  they 
can  be  produced  from  the  PI  method. 

4.  The  work  previously  described  considered  Iterative  procedures 
to  satisfy  various  constraints.  More  effort  Is  required  to 
handle  generalized  situations  involving  prescribed  volume, 
surface  area,  and  multiple  design  points.  From  this  stand¬ 
point,  viscous  effects  for  example  such  as  those  dealing 
with  friction  drag,  adverse  pressure  gradients,  are  all 
areas  that  are  of  interest. 

5.  For  both  the  steady  and  unsteady  problems,  the  materialization 
of  shocks  In  the  near  field  as  well  as  special  asymptotic 
boundary  layers  required  to  treat  these  structures  are  of 
Interest  regarding  matching  In  higher  orders  of  the  slender 
body  theory.  These  aspects  are  Important  with  respect  to 
nonlinear  lift  corrections  in  the  slender  body  theory,  and 
studies  of  the  Impact  of  lift  on  the  equivalence  rule  as 
described  in  Refs.  11  and  12. 
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PART  2 

ASYMPTOTIC  THEORY  OF  SOLID  TUNNEL  WALL  INTERFERENCE 
ON  TRANSONIC  AIRFOILS* 

SUMMARY 

The  method  of  matched  asymptotic  expansions  is  utilized  to  study 
the  singular  perturbation  problem  of  solid  wall  Interference  on  transonic 
airfoils.  For  moderate  to  large  wall  heights,  the  (Inner)  near  field  Is 
represented  as  a  linear  perturbation  about  the  nonlinear  free  field  which 
Is  assumed  to  be  governed  by  the  Karman-Guderley  small  disturbance  theory 
which  Is  non-unlformly  valid  as  the  walls  are  approached.  In  the  far 
field  (outer)  region,  another  approximate  representation  of  the  wall- 
airfoil  Interaction  Involving  a  multipole,  dominated  by  a  vortex  reflected 
between  the  walls  Is  valid.  Through  the  use  of  intermediate  limits, 
matching  of  both  representations  Is  demonstrated.  Some  numerical  solu¬ 
tions  for  the  Inner  problem  ?*>e  Illustrated  in  which  the  inner  limit  of 
the  outer  solution  Is  employed  as  a  far  field  boundary  condition  for  the 
perturbed  flow.  Means  of  correcting  the  tunnel  Incidence  to  obtain  an 
Interference- free  value  for  the  lift  are  demonstrated  from  the  examples. 

By  virtue  of  the  nature  of  the  perturbation  method,  the  height  dependence 
Is  separated  out  from  the  problem  and  universal  correction  functions  are 
available  from  the  theory  for  airfoils  at  given  Incidence  and  Mach 
number  conditions. 


*To  be  presented  as  AIAA  paper  82-0933  at  the  AIAA/ASME  3rd  Joint  Thermo¬ 
physics,  Fluids,  Plasma  and  Heat  Transfer  Conference  June  7-11*  1982/ 

St.  Louis,  Missouri. 
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1.0  INTRODUCTION 

Obtaining  free-fleld  characteristics  of  transonic  wind  tunnel 
models  will  continue  to  be  of  central  Importance  to  aeronautical  technology 
for  the  indefinite  future.  Free  flight  aircraft  performance  predictions 
depend  on  this  process  which  Is  complicated  by  the  presence  of  complex 
nonlinear  phenomena  In  the  transonic  regime.  Although  the  wind  tunnel 
Interference  problem  has  received  considerable  attention  at  subsonic 
speeds,  and  a  number  of  classical  theories  have  been  developed,  e.g.. 

Refs.  1-3,  there  Is  room  for  progress  even  In  this  linear  regime,  as 
exemplified  by  the  recent  efforts  of  Kraft,"  who  devised  an  Ingenious 
procedure  based  on  Cauchy's  Integral  formula  for  two-dimensional  flow 
over  confined  airfoils.  Among  the  advantages  of  this  method  In  contrast 
to  the  older  concepts  Is  that  by  measurement  of  two  flow  variables  It 
eliminates  the  need  for  both  knowledge  of  the  wall  characteristics  and 
analytical  synthesis  of  the  model.  Presently,  effort  Is  underway  to 
extend  this  procedure  to  three-dimensional  and  nonlinear  transonic  flows. 

In  connection  with  the  transonic  regime,  other  procedures  have 
been  developed  which  possess  attractive  features.  One,  which  Is  of  great 
Interest,  Is  a  scheme  that  has  been  developed  by  Murman.  This  "post  test 
assessment  method",  which  has  been  Implemented  In  two-dimensions  in  Ref.  5, 
Is  similar  In  some  respects  to  an  approach  developed  earlier  by  Kemp,6 
and  assumes  a  knowledge  of  the  experimental  pressures  on  the  model  and 
walls.  It  uses  modern  computational  and  optimization  procedures  to  deter¬ 
mine  whether  the  tunnel  Mach  number  and  model  angle  of  attack  are  correc¬ 
table  In  the  sense  that  almost  free-fleld  model  pressures  can  be  obtained 
at  practically  altered  values  of  these  parameters.  For  this  purpose,  an 
Inverse  problem  needs  to  be  solved.  In  three-dimensions,  surface  pressures 
over  the  model  are  generally  not  available,  and  efforts  are  currently  under 
way  to  modify  the  method  toward  the  use  of  less  Information  regarding  the 
model  near-field.7  If  this  goal  can  be  achieved,  this  process  will  be  of 
value.  Including  treatment  of  cases  in  which  the  walls  are  relatively 
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close  to  the  model.  In  spite  of  the  potential  utility  of  this  method, 
there  Is  a  need  for  approaches  that  can  reduce  the  number  of  Input 
parameters  necessary  to  compute  the  correction,  shed  light  on  the  physics 
of  the  wall  Interference  phenomena,  simplify  the  necessary  computations, 
and  be  generalized  to  three-dimensions,  as  well  as  unsteady  flows.  Asym- 
totlc  procedures  provide  such  advantages.  Furthermore,  they  can  provide 
valuable  Interactions  with  the  other  methods  previously  mentioned  to 
suggest  possible  Improvements  as  well  as  deriving  beneficial  features 
from  them  as  well. 

In  this  report,  an  application  of  the  asymptotic  method  to  the 
two-dimensional  case  Involving  a  transonic  airfoil  between  solid  walls 
will  be  studied.  This  effort  is  a  forerunner  of  three-dimensional 
formulations  which  besides  their  own  intrinsic  utility  have  relevance  to 
the  methods  of  Ref.  4  and  5.  Regarding  the  latter,  one  concept  presently 
being  considered  to  simplify  its  three-dimensional  application  is  unfold¬ 
ing  the  singular  character  of  the  near-field  rather  than  obtaining  the 
model's  shape  from  an  Inverse  solution  of  a  problem  Involving  specified 
pressures  based  on  measured  values.  As  will  be  evident  from  what  follows, 
the  asymptotic  procedures  can  provide  this  singular  behavior  from  certain 
limit  processes.  Moreover,  nonlinear  Integral  theorems  as  well  as  the 
asymptotic  structure  of  nonlinear  Integral  equations  arising  In  the  match¬ 
ing  scheme  occurring  In  the  asymptotic  analysis  could  be  of  use  In  the 
procedure  of  Ref.  4. 

Whereas  methods  such  as  those  of  Refs.  1-7,  can  handle  arbitrary 
model  to  wall  dimension  ratios  e,  the  perturbation  procedure  assumes  e  to 
be  small.  This  approximation  Is  useful  for  many  practical  cases.  Further¬ 
more,  even  for  situations  where  e  Is  not  so  small,  the  expansions  appear 
to  have  extended  validity. 

A  previous  analysis  along  these  lines  was  conducted  by  Chan' 
who  treated  the  two-dimensional  porous  wall  case.  Because  of  Its  Interest 
In  connection  with  compliant  boundary  applications  and  the  fact  that  the 
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Chan  solutions  apparently  do  not  subsume  It  as  might  be  expected  by  a 
taking  a  limit  of  the  porosity  parameter,  we  analyze  the  solid  wall 
case  In  this  report.  By  contrast  to  the  method  of  Ref.  8,  we  employ 
"Intermediate  variables"  to  match  the  model  near  and  far  fields. 

Although  slightly  more  cumbersome  than  the  approach  utilized  by  Van  Dyke,9 
we  believe  It  provides  a  reliable  means  of  ensuring  that  all  the  proper 
terms  are  Included  in  both  representations. 

In  addition,  numerical  results  will  be  provided  In  the  report  to 
give  some  Insight  into  the  nature  and  magnitude  of  the  interference  effects. 
The  work  of  Ref.  8  in  this  sense  was  strictly  formulation,  with  no  compu¬ 
tational  application  or  quantification  of  the  interference  given. 
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2.0  ANALYSIS 

Referring  to  Fig.  1,  which  Indicates  an  airfoil  nearly  midway 
between  solid  walls.  In  which  £  -  y6w,  x  and  y  are  Cartesian  coordinates 
In  units  of  the  chord,  h  Is  the  tunnel  half  height  In  the  same  units,  UT 
and  Mj  are  respectively  free-stream  speed  and  Mach  number  In  the  tunnel, 
and  transonic  flow  conditions  are  assumed  In  accord  with  a  Karman  Guderely 
(KG)  limit.  Thus,  If  6  Is  a  characteristic  thickness  ratio  of  the  airfoil 
whose  shape  Is  given  by 

*  *  iFg,*<x)  *  V-  1Fu.«Mflxl  '  1  (1> 


where  the  thickness  ratio  Is  considered  as  unvarying  from  the  (unbounded) 
free  field*  and  subscripts  u  and  it  denotes  upper  and  lower  surfaces  and 
ay  Is  the  angle  of  attack  In  the  tunnel,  the  appropriate  expansion  for 
the  velocity  potential  Is 

$  =  Uj  Jx+(5 2/3  $(x,y;  Kj,  Ay,H)  +  ...j  (2a) 


which  Is  valid  In  the  KG  limit 


1-Hi  aT  1/3 

x,y,  K.  «  — —  ,  A  ■  -p,  H  »  h6  fixed  as  6-*0, 

I  62/3  6 


(2b) 


*More  general  situations  could  be  considered. 
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and  leads  to  the  following  small  disturbance  problem*  for  the  confined 
transonic  airfoil: 


[KT  ’  (y+1>*x]  *xx  +  *  0 

(3a) 

$j(x,Q±) 

■  Fuy(*)  - 

(3b) 

<J^(x,H±) 

*  0 

(3c) 

r*j 

«  0,  Kutta-Joukowskl 

(3d) 

x  TE 
(x»l) 

Condition 

2.1  Inner  Problem 

A  simplification  of  the  problem  associated  with  Eqs.  (3)  can  be 
considered  in  which  expansions  of  $  In  (2a)  are  considered  for  H  -►  ». 

Near  the  airfoil,  the  walls  are  remote  and  thereby  Introduce  weak  pertur¬ 
bations  about  a  dominant  nonlinear  (KG)  free  field.  Accordingly,  an 
Inner  expansion  is  valid  which  can  be  written  as 


♦ 


$0(x,y)  +  ~  ^(x,*)  +  ... 
H 


(4a) 


x,J  fixed  as  (H  *  e"1)  ♦  »  (Inner  limit)  (4b) 


In  which  subscript  0  denotes  free  flight  and  1  is  the  dominant  wall 
correction.  The  gauge  function  ^  Is  anticipated  as  shown  to  achieve 

*In  the  Results  section  this  will  be  termed  the  "exact"  problem. 
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asymptotic  matching  to  be  discussed  subsequently.  For  generality  and 
use  In  tunnel  assessment  and  correction  procedures*  the  similarity 
parameter  and  angle  of  attack  are  also  written  as  perturbations  about 
the  free  field  values  along  with  those  for  the  potential  In  (4).  Thus, 


KI  *  KF  ♦  "c 


AT  *  AF  +  H2'  AC  + 


(5a) 


(5b) 


where  subscripts  F  denote  the  free  field,  C  is  associated  with  wall 
corrections,  and  T  designates  the  tunnel  value.  On  substitution  of 
(4)  and  (5)  into  the  problem  (3),  and  retaining  corresponding  orders, 
the  zeroth  and  first  order  problems  are  obtained  as 


IKp  -  (y+1)$q  j  $q  +  =  0 

'  x/  xx  yy 


l^p  "  (Y+1)$Q  "  (Y+l)it!i  ^  $1  -  “Kg  <f> 

'  X'  XX  *  **  Av/\/ 


4\{x,0 t)  *  F'Ut£(x)  -AF 

4>q  (x,0±)  •  -Ac 


WnWn- 


x  xx  yy 


X  < 


xx 
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In  the  problem  associated  with  the  "variational"  equation  (6b)  which 
resembles  others  such  as  the  lifting  line  theory  of  Ref.  10,  the  jumps 
across  the  0^  order  shocks  are  perturbed  according  to  the  developments 
given  therein.  In  (6),  if  r  Is  the  circulation,  it  can  also  be  written 
in  the  perturbation  form 


To  complete  the  formulation,  a  far  field  Is  required.  In  this  connection, 
it  is  crucial  to  note  that  the  approximation  (4a)  is  non-uniformly  valid 
near  the  walls  and  an  "outer"  representation  is  required.  An  "inner" 
limit  of  this  outer  solution  (approaching  the  boundaries  of  the  airfoil 
near  field)  supplies  the  required  inner  far  field  boundary  conditions. 
These  will  be  discussed  presently.  For  the  far  field  of  the  0th  order 
problem  the  representation. 


*0 


~rF9  /  rF  \2y+l  £n  r 
2tt  V  2tt  /  4Kp  r 


cose 


/  (v+1)  cos 36 

\2tt/  16Kp  r 


+  -fipsfi  +  fo  ..?j,n,e 
2tt  r  2tt  r 


(7a) 


*7x2+  Kpy: 


,  jCy 

tan*1  ~ —  fixed,  (ray  limit)  (7b) 


Is  appropriate  on  the  basis  of  Ref.  10.  In  addition,  the  far  field 
anticipated  and  to  be  refined  for  is 
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BjXy  +  CjX 


V 


♦  N, 


rc\ 

TT 


(7c) 


In  which  x2,^2  quadratic  terms  are  missing  based  on  the  matching  to 
be  Indicated,  and  and  Nj  are  constants  to  be  determined  In  this 

matching. 


In  this  simulation,  the  aim  Is  to  select  Ac 
the  correction  In  some  sense. 


and  Kg.  to  minimize 


2.2  Outer  Problem 

At  distances  from  the  airfoil  where  the  wall  Influence  Is  not 
small,  the  outer  expansions 


♦  ■  *0(x*,y*)  +  ^  ^(x*,y*)  +  £  ^(x*,y*)  +  ... 


(8a) 


are  anticipated  as  valid  In  the  limit 

x*  *  y*«  ^  fixed  asH-»»,  (outer  limit),  (8b) 


where  the "switchback"  term  Is  Inserted,  anticipating  In  the  matching 
the  effect  of  logs  arising  from  the  outer  limit  of  (7a).  On  substitution 
In  (3),  the  equations  obtained  for  the  terms  in  (8a)  are 

Mb  +  *  0  ,  0(1/H2)  (9a) 
f  Vx*  yy*y* 
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r  *X*X*  iy*y* 


KF^1  + 

lx*x*  Ay*y* 


0  ,  o(*nH/HJ) 

(y+1)*q  *>q  » j(i/h*) 

ux*  ux*x* 
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(9b) 

(9c) 


with  the  wall  boundary  conditions. 


(x*,±l)  *  (x*,±l)  *  (x*,±l) *  0.  (9d) 

uy*  *y*  Ay* 

2.3  Preliminary  Matching  Considerations 

Under  the  assumption  that  to  dominant  order,  V’g  is  a  potential 
vortex  reflected  between  parallel  walls,  the  exact  expression  from  a 
Prandtl-Glauert  scaled  version  of  incompressible  flow  theory  Involving 
either  images  or  conformal  mapping  is 

+  1ip0  *  F(Z)  =  complex  potential 

*  in  tanh  (10a) 

Z  -  -4l-  +  1y*  =  X  +  iY  (10b) 

vK  f 


The  validity  of  Eqs.  (10)  Is  predicted  on  the  assumption  that  the  far  field 
is  subsonic.  On  expansion  of  (10)  for  Z  ■+  0,  it  is  evident  that 
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\  •  £  ]-W  ^  «V  ♦  ...  J  .x‘OT-0  .  (11)  (11) 

In  (11),  the  x*y*  term  Is  the  effect  Induced  In  the  origin  due  to  the 
wall  reflections. 

To  compute  symmetry  properties  that  can  be  Inferred  from 
the  right  hand  side  of  (9c)  are  useful.  These  can  be  obtained  from 
symmetries  of  which  are 

<F0(x*,y*)  *  -*0(x*,-y*) 

*0(x*,y*)  *  -*Q(-x*,y*) 

from  which  it  follows  that 

^(x*,y*)  ■  ^(x*,-y*) 

^(x^.y*)  -  -fj(-x*,y*) 

With  these  relations  In  mind,  the  slmular  behavior  required  for  matching 
can  be  obtained  from  a  linear  combination  of  particular,  complementary, 
and  homogeneous  solutions  of  (9c),  In  which  the  complementary  solution 
Is  added  to  satisfy  the  boundary  conditions,  and  the  homogeneous  solu¬ 
tions  are  required  to  provide  enough  arbitrariness  In  the  doublet 
strengths  to  achieve  matching.  Denoting  as  a  particular  solution, 
a  choice  Is  made  for  <pp  In  which 


%  ^0  * 
°X  °XX  *F  uY  °XY 


(12) 
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where  the  right  hand  member  of  (12)  Is  obtained  from  the  Cauchy-Rlemann 
equations  and  Is  used  In  favor  of  the  middle  one  to  avoid  branch  point 
singularities  at  the  origin  such  as  that  occurring  with  the  function 
<8>-  arg  Z,  and  4»0  Is  the  harmonic  conjugate  of  Using  a  clever 
device  employed  In  Ref.  11,  the  appropriate  particular  solution  of  (12) 
can  be  verified  to  be 


*P "  { Vox  *  Y*ox\} 


On  expanding  (13)  near  the  origin  Z*0,  some  of  the  necessary  terms 
required  for  matching  with  (7a, c)  appear,  l.e.. 


1>  “  cos9  *  \  *  [£n(?7^,)+ 

•••}  (15a) 


as  r*  ■  Vx^+K^y**  =  r/H  -*■  0 


(15b) 


e*  -  tan  -£*— 


tan  — —  *  e,  fixed. 


(15c) 


The  singular  behavior  of  (14)  near  the  origin  contains  no  branches  as 
desired.  Consistent  with  earlier  remarks,  the  appropriate  doublet  terms 
Involving  are  apparent  and  are  of  specific  strength.  However, 
additional  ones  of  arbitrary  adjustable  strength  are  needed  for  matching. 
These  will  be  obtained  from  a  suitable  homogeneous  solution  of  (9c). 
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In  view  of  these  considerations,  the  complete  solution  of  (9c)  Is  written 
as 


*1 


V 


(16) 


where  V?  Is  a  harmonic  complementary  solution  added  to  satisfy  the 
boundary  conditions  violated  by  According  to  (9d), 


3*C(X.*1)  3*„(X,±1)  _  (r*l)  /rF\2  tt3  shirX 

ir  mmw  *  r  (Rhixjf 


(17) 


A  Fourier  integral  representation  for  v>  can  be  obtained  from  the  wall 

-3*(X,±1)  c 

source  distribution - -  in  (17)  as 


*c  %  i  (^)  r  J[ dk  £cosk(x’^  d+Sfe  ^ 


which  upon  Integration  simplifies  to 


-  _  (y+1)  /  rF\zl  f* sinkX  chkY  ... 

*c  *  JJ,  iTT^ —  dk 


The  behavior  of  which  as  X,  Y  0  is 


(18) 


To  obtain  In  (16),  a  further  decomposition  In  which 
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Is  used,  where  ^  and  are  arbitrary  strength  doublets  that  are 
respectively  symmetric  and  antisymmetric  In  y*  and  <*>sc  and  are 
complementary  solutions  added  to  each  of  these  components  to  satisfy 
the  wall  boundary  conditions.  Thus, 


d0  cose* 
Er"7*"“ 


d0  x* 
Tn 


(20a) 


*0  sine* 
F  r* 


e0  v* 

F  xfi+Kpy*5  ’ 


(20b) 


with  dQ  and  eQ  the  adjustable  strengths  to  be  determined  from  matching. 
Furthermore,  on  using  the  Fourier  Integral,  recognizing  that 

3<&r(X,±l)  dn 

TT  *I2^T1W 

3^AC(x.iD  eQ  1-X2 

ir  “Fv^UW  * 

the  following  representation  of  the  complementary  components  ^  Is 
obtained: 


(20c) 
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In  a  similar  manner. 


V'i 


£dk  £C«k'X-'>  {te)Sir^fi*5=l!}d£ 


“JnJtjJI  **  “  r*  *° 


(2W) 


where  the  antisymmetric  complementary  component,  represents  a 
dl vortex  at  the  origin,  which  gives  an  Induced  upwash. 

2.4  Matching 

With  the  presumption  that  there  Is  an  overlap  domain  of 
validity  of  the  inner  (near  field)  representation  of  $  and  Its  outer 
(far  field)  form,  with  no  intermediate  expansion  required,  appropriate 
developments  in  this  domain  for  each  are  to  be  obtained  In  what  follows. 
From  these  cornnon  representations,  the  unkown  elements  in  both  are 
determined  by  matching.  Coordinates  specifying  the  overlap  domain  are 
obtained  from  the  ray  preserving  intermediate  limit 

\mnm  •  ’  dir  •  m‘d  “  H  *  *  • 

where  n(H)  ■*  »,  is  an  order  class  Intermediate  between  unity  and  H. 

Thus,  x*  ■  x^  -*■  0,  y*  ■  jj  y^  -*■  0,  and  x  *  nx^  ■*•«*.  £  ■  ny^  ■*  •. 
nr 

Also  r*  ■  with  r*  ■  x*  +  K^y*.  The  matching  Is  effected  by  writing 
the  Inner  and  outer  representations  in  terms  of  these  Intermediate 
variables.  For  the  outer  expansion,  this  Is 


* 
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♦  *  +  +  •••  +  TT*i(?  Syr) 


♦A 


nr 

*°-r 


0 


cose 


0 

1  V±1  /  rp  V  J  1  COS  36 

IT 


0 

,2 


TT 

(8)  0 


w*  nxn  *2  n\ 


■  if ln  TJXfwi  *  rsjfl cos2e  +  ••• 

(8)  0 


1  &i(iViiISu  +1^0 

i^Wi2ir +  —  +Tfi? 


cose 

nr 


0 


@ 


1  d0  TT2  nX~  .  1  eQ  sine  +  1  ®0  Tr2  +  (n) 

*JT5XjT2'T*  hUT—  hTQ^^-T  *  •••  *  <21' 

nr 

where  similar  terms  have  been  identified  with  identical  numerical 
labels  for  matching  purposes.  From  the  orders  appearing  In  (21)  a 
sharper  lower  bound  for  n(H)  is  obtained  giving 


«  n  «  H 


For  the  inner  problem,  If  a  new  dependent  variable  Is  Introduced, 
in  which 
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,  %  *c 

5  yff  *  +  *1  • 

then  (6b)  gives 

(My+1)<|>0  )$,  -  (y+1)$,  $n  +  *.  »  0  , 

x  Axx  Ax  yxx 

and  on  the  basis  of  the  foregoing  matching  requires  that 

rCe 

$j(x,y)  *  B^xy  +  P^xinr  ♦  Cj  x  ♦  Q^x  cos28  ♦  +  ... 


or 


^(x,?)  •  -BjXjf  +  PjXtnr  +  +  C^x+QjX  cos26 

f* 

♦V  +  ni-tS,+  —  (22) 

In  the  far  field  r  ■  Vxa+Kpy*  -*■  •,  In  which  the  constants  Bi*pl*cl* 
HjtHj  and  Qj  are  to  be  determined  by  the  matching  process.  For  this 
purpose,  and  Inner  solution  Is  written  In  Intermediate  variables.  The 
resulting  expansion  is 


0 

+  °°  cose 

7T  nr„ 
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+  *0  sine 
IFF  nrn 


+  °Hr)*  W  j  BlnVn  *  P1TV"  r,rn  +  (wH  +  Cl)  nxn 
<S>  © 

rc  ) 

♦  Qjnx^  cos29  ♦  Mjny^  -  ^  e  +  ...  |  .  (23) 


In  (23),  the  terms  labeled  Q),(2),and  (3)  match  automatically  to  their 
similarly  labeled  counterparts  In  (21)  If  the  harmonic  switchback  term 
is  given  by 


„  ,  y+1  M  COse 

%  w  ~  • 


On  the  basis  of  comparison  of  the  orders  to  be  indicated,  the  following 
matchings  apply: 


(ip) :  8i  ■  w (s?) ^7  •  ®  -(&)  '■  pi  ’  (sr) (Mb) 


■  ■  ••  v  •  •  . 
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(24c) 


(2)  =»  qi  *  - 

(24d) 

©  =>  M1  *  JQTf  IS 

(24e) 

To  complete  the  matching  considerations,  it  is  noted  that  a  switchback 

term  is  required  to  match  the  term  $2  2§^~  x^fcnH  arising  from 

(Z>  in  the  outer  expansion.  Since  this  £wi tchback^term  does  not 
vanish  at  infinity,  a  corresponding  switchback  term  can  be  added  to  the 
inner  solution  with  a  corresponding  shift  in  the  transonic  similarity 
parameter.  Accordingly,  the  revised  inner  expansions  replacing  (4a) 
and  (5a)  are 


*  $0  ♦  "pr  +  W?  +  •  •  •  (4a‘) 

Kf  *  Kp  +  +  ft?  +  •  •  •  •  (5a') 

To  match,  therefore,  Is  selected  as 
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which  satisfies  the  following  switchback  equation  obtained  by  substitution 
of  (4a*)  and  (5a')  Into  (3a): 

<Kf'tY+1)V\,  * (v+1>  V°*x +  \  *  ■**v  • 

providing 


K 


S 


(26) 


This  completes  the  matching.  All  the  constants  are  given  Eqs.  (24),  (25) 
and  (26)  with  the  exception  of  the  doublet  strengths  Dq  and  Eq.  These 
are  obtained  from  an  extension  of  the  Integral  equation  analysis  given  In 
Ref.  12.  Although  the  dominant  term  Is  given  in  Ref.  13,  an  analysis  Is 
required  to  obtain  the  doublet  contributions  which  are  higher  order.  The 
details  are  contained  In  the  Appendix  and  the  final  results  are 


(27a) 

(27b) 


\ 5  J £<  w* 

with  p  and  6  as  polar  coordinates,  and  £  Is  a  dummy  variable  for  x. 
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3.0  RESULTS 

To  obtain  the  wall  interference  corrections  to  be  discussed 
in  what  follows,  the  free  field  problem  (6a,c,e)  was  solved  by  the 
SLOR  procedure  of  Ref.  12.  The  perturbation  problem  associated  with 
(6b,d,e)  with  the  far  field  (22)  was  solved  with  the  linear  counter¬ 
part  of  this  method  using  modifications  of  special  procedures  to 
handle  the  perturbed  shock  waves  developed  by  R.D.  Small11*  to  compute 
solutions  of  the  lifting  line  problem  formulated  in  Ref.  10. 

As  an  indication  of  the  size  of  the  basic  effect,  surface 
pressures  over  a  NACA  0012  airfoil  at  Mp  =  0.7,  ap  =  3.5°,  were  computed 
for  A  =  Kc  *  0  for  H  *  1,  or  h*2.  The  perturbation  values  are  com¬ 
pared  against  those  of  the  free  field  in  Fig.  2  and  show  a  surprisingly 
small  interference  effect  for  this  small  H  value.  For  this  and  larger 
H,  a  part  of  this  trend  is  associated  with  the  small  coefficients  of  the 
gauge  functions  of  H  appearing  in  (4a).  A  relevant  factor  may  be  the 
inability  of  the  theory  to  generate  a  shock  from  an  initially  shockless 
flow.  As  in  transonic  lifting  line  theory,  only  shockless  perturbations 
about  shockless  flows  or  shocked  alterations  about  shocked  cases  are 
accessible  to  this  theory.  For  the  latter,  the  jump  conditions  given 
in  Ref.  10  must  be  used  for  the  linear  variational  problem  ^ ,  as  pre¬ 
viously  indicated.  When  shocked  perturbations  occur  on  shockless  flows, 
a  special  inner  solution  near  the  shock  Is  required.  In  spite  of  these 
limitations,  it  is  felt  the  theory  of  this  report  could  provide  useful 
information  for  the  typical  case  in  which  shocked  flows  remained  shocked. 

In  Fig.  3,  free  and  confined  chordwise  pressures  over  the  same 
NACA  0012  airfoil  for  different  heights  at  a  higher  Mach  number, 

Mp  *  0.75  are  shown.  At  the  incidence  indicated  (op  ■  2°),  the  typical 
supercritical  pressure  distributions  are  evident.  As  anticipated,  the 
confined  pressures  approach  the  free  field  values  as  H  Increases. 
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x 


Fig.  2.  Free  and  confined  chordwlse  pressures  NACA  0012  airfoil 
Mc"0.7,  ac-3.5°,A  -K  -0,  H»1 
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.  Free  and  confined  chordwl se  pressures  HACA  0012  airfoil 
Hf-0.75,  op-20,  Ac-Kc-0 


2-24 


Rockwell  International 

Sc  tone*  Center 


SC5267.3FR 

A  comparison  of  the  perturbation  solution  of  this  report  with 
that  of  the  "exact"  problem  embodied  in  Eqs.  (3)  Is  shown  In  Fig.  4. 

The  agreement  is  only  fair  for  this  low  value  of  H  with  considerable 
degradation  downstream  of  the  shock.  This  may  be  related  In  some  degree 
to  the  numerical  approximations  required  to  handle  nodes  in  the  vicinity 
of  the  shock,  and  the  need  for  a  local  solution  there. 

As  anticipated,  the  agreement  between  exact  and  perturbation 
solutions  improves  with  an  increase  in  K  as  shown  in  Fig.  5.  However, 
here,  there  appears  to  be  some  deterioration  near  the  trailing  edge.  In 
making  these  comparisons,  a  certain  degree  of  sensitivity  and  nonunique¬ 
ness  was  experienced  with  respect  to  restarts  of  the  exact  and  perturbation 
codes  to  obtain  the  various  parametric  ases.  Multiple  solutions  and 
hystereMc  behavior  have  been  ascribed  to  underlying  properties  of  the 
continue,,  problem  in  the  free  field  case  for  the  full  potential  equation 
in  Ref.  15.  The  confined  problem  may  also  have  these  features.  Since 
small  differences  are  involved  in  the  interference  assessments,  under¬ 
standing  the  source  of  these  nonuniquenesses  deserves  attention  in 
future  effort  not  only  with  asymptotic  methods  but  in  purely  computational 
tunnel  wall  procedures  as  well.  Other  aspects  associated  with  validations 
of  this  type  relate  to  truncation  errors  of  the  computational  solutions, 
the  role  of  the  far  field,  initial  conditions  for  the  SLOR  algorithm  and 
the  degree  to  which  the  Kutta  condition  is  satisfied. 

As  a  further  indication  of  the  significance  of  some  of  these 
factors,  the  lift  coefficients  for  the  perturbation  and  exact  solutions 
have  been  compared  at  various  heights  for  the  two  cases  considered  in 
Figs.  1-4.  In  Table  1,  excellent  agreement  is  shown  over  the  entire  range 
of  heights  with  some  understandable  increased  discrepancy  at  H  *  1.  In 
view  of  the  accuracy  factors  previously  enumerated,  this  agreement  is 
somewhat  better  than  would  be  anticipated. 
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.  Comparison  between  exact  and  approximate  chordwlse  pressures  on 

confined  airfoil  NACA  0012  airfoil  Hc»0.75,  V*^0,  H«l, 
a  -1/  «  r 
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Comparison  between  exact  and  approximate  chordwl se  pressures  on 
confined  airfoil  NACA  0012  airfoil  HF«0.75,  ap»2°,  H-2, 
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TABLE  1.  COMPARISON  OF  LIFT  COEFFICIENTS  C.  FOR  CONFINED 
NACA  0012  AIRFOIL  AT  VARIOUS  WALLlHEIGHTS, 

Mf  -  .7,  otp  *  3.5° 


H 

Methoa^^ 

1 

2 

3 

4 

5 

Exact 

Asymptotic 

Solution 

m 

■ 

m 

m 

m 

In  Table  2,  a  similar  comparison  is  shown,  again  indicating  good 
agreement,  but  with  increasing  discrepancies  at  the  lower  H's. 


TABLE  2.  COMPARISON  OF  LIFT  COEFFICIENTS  C.  FOR  CONFINED 
NACA  0012  AIRFOIL  AT  VARIOUS  WALLlHEIGHTS, 

Mf  =  .75,  dp  =  2° 


H 

MethoiT\^^ 

1 

2 

3 

4 

5 

6 

Exact 

Asymptotic 

Solution 

0.4322 

0.4356 

0.4256 

0.4264 

0.4234 

0.4248 

0.4234 

0.4242 

0.4233 

0.4238 

0.4234 

0.4237 

It  is  interesting  to  note  that  strictly  In  accord  with  the 
asymptotic  order  of  magnitude  of  the  second  order  switchback  term  of 
Eq.  (4a'),  the  perturbations  associated  with  the  solid  walls  are 
0(tnH/H2)  which  are  smaller thas  the  0(H“l)  one  for  the  perturbation 
solution  relevant  to  the  porous  case  described  In  Ref.  8.  Since 
porosity  Is  Introduced  to  weaken  the  wall  Interactions,  this  would 
superficially  seem  paradoxical.  However,  the  coefficients  of  the  Indl 
cated  gauge  functions  of  H  do  also  play  a  role.  From  a  numerical 
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viewpoint,  these  could  effect  the  trend.  Another  perhaps  more  salient 
consideration  deals  with  one  motivation  for  porosity  as  a  means  of 
cancelling  compression  and  expansion  wave  reflections.  This  mechanism 
Is  of  limited  relevance  to  the  subsonic  far  field  assumed  In  the  formula¬ 
tion  of  this  report.  Accordingly,  It  Is  reasonable  to  state  that  whereas 
perforated  walls  can  minimize  supersonic  wave  reflection  Interference, 
they  may  not  be  optimum  to  reduce  the  subsonic  constriction  type.  This 
aspect  Is  probably  related  to  the  order  of  magnitude  considerations 
mentioned  previously. 

Regarding  Interference  minimization.  Fig.  6  shows  the  interference 
lift  as  a  function  of  the  tunnel  wall  Interference  angle  of  attack 
parameter  A~.  The  Intersection  with  the  absissca  Is  the  tunnel  Incidence 
correction  required  to  produce  the  free  field  lift  characteristics.  Fig.  7 
gives  a  similar  curve  at  a  larger  Mach  number  but  smaller  geometric  Inci¬ 
dence  requiring  a  larger  Ac  correction.  In  accord  with  (4a'),  variations 
of  this  type  can  be  used  as  a  universal  curve  valid  for  arbitrary  H's 
within  assumptions  of  the  asymptotic  theory.  A  similar  conclusion  applies 
to  the  other  forces  and  moments  as  well  as  the  surface  pressures.  For 
pressures,  the  Interference  can  only  be  reduced  In  a  mean  squared  sense 
by  adjusting  Ac  or  Kj.. 
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4.0  CONCLUSIONS 

An  asymptotic  theory  of  solid  tunnel  wall  Interference  on 
transonic  airfoils  has  been  developed  based  on  the  limit  of  large  wall 
height  within  the  Karman  Guderley  (KG)  theory.  This  approximation 
leads  to  a  singular  perturbation  problem  In  which  the  near  field  approxi¬ 
mation  of  weak  linear  perturbations  about  the  KG  small  disturbance  model 
Is  non-unlformly  valid  as  the  walls  are  approached.  A  far  field  for 
the  near  field  problem  comes  from  an  "outer"  approximation  Involving 
wall  reflections  of  a  multipole  representation  of  the  airfoil  dominated 
by  a  reflected  potential  vortex.  From  the  analysis  and  results,  the 
following  observations  can  be  made: 

e  Switchback  terms  are  needed  In  the  transonic  case  to 
properly  match  the  near  and  far  fields.  These  alter  the 
effective  similarity  parameter  of  the  confined  flow  field. 

•  Perturbations  of  the  pressures  and  forces  in  the  near  field 
are  0(£nH/H2). 

e  The  asymptotic  theory  provides  good  agreement  with  the 
exact  KG  model  for  large  to  moderate  H. 

e  The  usefulness  of  the  theory  is  enhanced  by  the  fact  that 
It  separates  H  out  of  the  Interference  problem  and  requires 
the  solution  of  a  nonlinear  and  linear  problem  rather  than 
two  nonlinear  problems  representing  respectively  the  free 
and  confined  fields  to  obtain  the  Interference  correction. 


The  procedures  developed  In  this  report  provide  a  basis  for  extension  to 
three-dimensional  and  unsteady  problems.  Lifting  line  and  slender  body 
theory  can  supply  the  necessary  starting  points  for  the  former. 
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APPENDIX 

Determination  of  Eq  and  Dq 

Using  Green's  theorem  In  a  manner  similar  to  that  given  in  Ref. 12, 
the  integral  equation  obtained  for  Is 

*<*•*>  ‘  45  '  rjf(^  45  *  1  •  (A1) 

where  the  subscripts  on  4>  have  been  dropped,  [f]w  *  fy=Q+  ~  fy=Q-*  a 
change  of  notation  from  the  previous  sections,  X  *  x,  Y  = 

S  *  i—  An  >/(x-E)*+(Y-n)5,  r  *  [<j>l  ,  E  and  n  are  dummy  versions  of 

fcTT  "«-1 


r 


* 


S^X.Y.E.n }$*(E.n)  d?dn 


(A 2) 


Ref.  12  is  restricted  to  the  nonlifting  case,  and  in  that  situation,  ^ 
dies  off  sufficiently  rapidly  so  t  tat  a  convergent  integral  involving 
the  average  source  strength  Is  obtained  by  placing  the  kernel  in  (A2) 
outside  the  integral  signs.  This  Is  true  since  $  in  the  nonlifting  case 
Is  asymptotically  a  doublet  and  Is  Oir'1)  as  r  ■*  ®.  By  contrast,  <|>  Is 
a  potential  vortex  in  accord  with  (7a)  and  can  no  longer  be  regarded 
as  approximately  of  compact  support.  To  deal  with  this  difficulty,  I 
is  written  in  the  following  form: 

1  "  !l  +  lz  +  h  *  *A3) 
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where 


I 


I 


S 


5 


(Rcose-pcosa) 
Rz-2Rocos(6-  4)+p'2 


X  ■  Rcose,  5  *  pcos*  . 


(A4) 

(A5) 


In  (A3),  I  has  been  regularized  at  «  by  subtracting  and  adding  the 
asymptotic  representation  of  4>^  in  that  neighborhood.  A  resulting  singu¬ 
larity  near  the  origin  is  avoided  by  decomposing  the  integral  into  a  part 
Inside  the  unit  circle  and  outside  of  it.  Since  its  integrand  is  of 
compact  support,  1^  can  be  approximated  by  the  method  of  Ref.  12,  l.e., 
approximation  of  the  integrand  for  large  R  and  integration  term  by  term. 
The  same  conclusion  applies  to  I2,  whose  Integrand  is  On^r)  as  R  ®. 

The  Integral  Ij  can  be  treated  exactly.  However,  the  operations  have 
to  be  performed  with  great  care  because  of  the  discontinuous  nature  of 
the  Integrand.  The  choice  of  the  lower  p  limit  as  unity  in  (A4)  and  (A5) 
leads  to  considerable  simplifications  as  compared  to  other  values  and 
therefore  was  adopted. 

The  results  for  Ij,  I2.  and  I3  are 


(A6) 
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‘a-twr^/o  <»> 

which  upon  substitution  in  (Al)  and  approximation  of  the  one-dimensional 
integrals  therein  leads  to  (27). 
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